GBE 



Patterning and Regulatory Associations of DNA Methylation 
Are Mirrored by Histone Modifications in Insects 

Brendan G. Hunt^ Karl M. Glastad^ Soojin V. Yi'^, and Michael A.D. Goodisman* 

School of Biology, Georgia Institute of Technology 

"Corresponding authors: E-mail: soojin.yi@biology.gatech.edu; michael.goodisman@biology.gatech.edu. 
■''These authors contributed equally to this work. 
Accepted: February 25, 2013 

Data deposition: Solenopsis invicta DNA methylation data have been deposited in Gene Expression Omnibus under the accession GSE39959. 

Abstract 

Eplgenetlc information is an important mediator of the relationship between genotype and phenotype in eukaryotic organisms. One 
of the most important and widely conserved forms of epigenetic information is the methylation of genes. However, the function of 
intragenic DNA methylation remains poorly understood. The goal of this study was to gain greater understanding of the nature of 
intragenic methylation by determining its role in the multilayered epigenetic landscape of insects. We first investigated the evolu- 
tionary lability of DNA methylation by examining whether methylation patterns were conserved in the fire ant and honey bee. We 
found that DNA methylation was targeted to largely overlapping sets of orthologs in both species. Next, we compared intragenic DNA 
methylation levels in the fire ant and honey bee to comprehensive epigenetic and gene-regulatory data from Drosophila melano- 
gaster orthologs. We observed striking evidence of a conserved association between DNA methylation in fire ants and honey bees, 
and several active histone modifications, constitutive gene expression, and "broad" promoter architecture in D. melanogaster. 
Overall, our study illustrates that DNA methylation is a single component of a conserved, integrated, multilayered epigenetic and 
regulatory landscape in insect genomes. 
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Introduction 

A wide array of animals and plants display intragenic DNA 
methylation. Intragenic methylation is the methylation of 
exons and introns (Zemach et al. 2010) and is targeted to 
active genes (Zemach et al. 2010; Jjingo et al. 2012), in con- 
trast to promoter methylation, which is linked to gene repres- 
sion (Bird and Wolffe 1 999; Pai et al. 201 1 ; Zeng et al. 201 2). 
Intragenic methylation has been hypothesized to affect a wide 
array of processes including the regulation of messenger RNA 
splicing (Shukia et al. 2011), the initiation of transcription 
within genes (Zilberman et al. 2007; Maunakea et al. 2010), 
and transcriptional elongation (Lorincz et al. 2004; Zilberman 
etal. 2007). 

Given the important regulatory roles suggested for intra- 
genic DNA methylation, it is somewhat perplexing that DNA 
methylation has been lost entirely in some eukaryotic lineages, 
including several clades of insects (Zemach et al. 201 0; Glastad 
et al. 201 1). Insight into the loss of DNA methylation may be 
found in the presence of functional overlap with histone 



modifications (Cedar and Bergman 2009; Nanty et al. 201 1; 
Suganuma and Workman 201 1). The regulatory roles of his- 
tone modifications include the mediation of binding affinities 
of protein complexes, such as those related to transcriptional 
and splicing machinery (Kolasinska-Zwierz et al. 2009; Luco 
et al. 2010, 2011; Negre et al. 2011). Thus, similar to DNA 
methylation, histone modifications convey important epige- 
netic information, which affects gene regulation and organis- 
mal development. Indeed, one of the most exciting 
developments in the field of epigenetics was the realization 
that patterns of DNA methylation interact with histone mod- 
ifications and nucleosome positioning in a multifaceted epige- 
netic landscape (Cedar and Bergman 2009; Chodavarapu 
etal. 2010). 

This study is focused on understanding the patterning and 
regulatory significance of intragenic DNA methylation in insect 
genomes. We first examine the extent to which DNA meth- 
ylation targeting is conserved between species. This provides 
insight into the evolutionary lability of DNA methylation. 
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Second, we determine whether conserved evolutionary 
associations exist between DNA methylation patterning and 
histone modifications. The associations between these epige- 
netic marks have remained largely unexplored in insects 
because the primary model for understanding insect histone 
modifications, Drosophila melanogaster, lacks DNA methyla- 
tion (Zemach et al. 2010). To address this issue, we compare 
DNA methylation data from the fire ant Solenopsis invicta and 
the honey bee Apis mellifera to extensive histone modification 
data from D. melanogaster. Finally, we explore associations 
between gene expression, DNA methylation, and histone 
modifications, and how these associations relate to other reg- 
ulatory features in insect genomes. Overall, our investigation 
provides important insight into how different layers of epige- 
netic information are linked to gene regulation. 

Materials and Methods 

DNA Methylomes 

We generated DNA methylation data from the fire ant 
5. invicta to compare with existing data from the honey bee 
A. mellifera (detailed materials and methods provided in 
Supplementary Material online). We used whole bodies of 
four different adult fire ant morphs (queens, workers, haploid 
males, and diploid males) to generate single-base resolution 
DNA methylation maps (DNA methylomes). From these data, 
we generated a species-level 5. invicta DNA methylome (sup- 
plementary table S1, Supplementary Material online). For 
comparative purposes, we mapped a species-level DNA 
methylome of A. nnellifera using adult queen and worker 
brain DNA methylomes (Lyko et al. 2010) and a whole body 
worker DNA methylome (Zemach et al. 201 0) (supplementary 
table S2, Supplementary Material online). 

Bisulfite conversion and lllumina paired end sequencing of 
5. invicta genomic DNA were performed by Beijing Genomics 
Institute (Zhenzhen, China). Reads from 5. invicta and A. nnel- 
lifera (Lyko et al. 2010; Zemach et al. 2010) were aligned to 
reference genomes (Honeybee Genome Sequencing 
Consortium 2006; Wurm et al. 2011) using Bismark 
(Krueger and Andrews 2011). Fractional methylation values 
were calculated for each CpG site as mCG/CGaii, where mCG 
is the number of reads with a methylated cytosine at a CpG 
site (according to nonconversion) and CGaii is the total number 
of reads mapped to the site. Fractional methylation values 
were then averaged across each annotated element. 
Solenopsis invicta DNA methylation data are available from 
Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/ 
geo/; GSE39959). 

Histone Modification Data 

Drosophila melanogaster preprocessed Chromatin Immuno- 
Precipitation using genomic tiling arrays (ChlP-chip) data files 
identifying regions of significant enrichment for each histone 



modification were obtained from modEncode (Celniker et al. 

2009) (supplementary table S3, Supplementary Material 
online). We used ChlP-chip data from the S2-DRSC late- 
stage embryonic cell line. We determined the proportion of 
an element overlapped by significant regions of enrichment 
(overlap), the mean enrichment score for these regions, and a 
composite metric, which is the product of score and overlap 
metrics. The composite metric was used in all analyses unless 
otherwise indicated. 

Gene Expression 

Solenopsis invicta gene expression data were obtained from a 
prior cDNA microarray study (Ometto et al. 201 1). For each 
gene, we assessed the coefficient of variation (standard 
deviation/mean; CV) of expression values as the mean of 
CV values calculated separately for whole body 5. invicta 
adult and pupal workers, queens, and haploid males 
(median of 5 biological replicates per morph) (Ometto et al. 
2011). 

Drosophila melanogaster gene expression data were ob- 
tained from FlyAtlas (Chintapalli et al. 2007). We estimated 
the overall CV in FlyAtlas expression values as the mean of CV 
values calculated separately for each of 1 0 tissues (based on 4 
biological replicates per tissue). The "tissue specificity index" 
of each gene was calculated using data from 10 tissues. 

Orthology 

We used OrthoDB (Waterhouse et al. 201 1) 12-insect orthol- 
ogy data to assign orthologous genes between 5. invicta, 
A. mellifera, and D. melanogaster. 

Gene Ontology 

Solenopsis invicta and A. mellifera Gene Ontology (GO) was 
assigned, and GO term enrichment determined, using 
Blast2G0 (Conesa et al. 2005). 

Results 

DNA Methylation Targets Are Conserved in the Fire Ant 
5. invicta and Honey Bee A. mellifera 

We found that DNA methylation was targeted primarily to 
exons in the fire ant and the honey bee (supplementary figs. 
S1-S3, Supplementary Material online) (Lyko et al. 2010; 
Bonasio et al. 2012). DNA methylation was lower in introns 
than exons, but the discrepancy between exon and intron 
methylation level was much larger in the honey bee (supple- 
mentary fig. S2, Supplementary Material online) (Lyko et al. 

2010) . We also observed a pronounced preference for DNA 
methylation targeting to 5^-intragenic regions in both species 
(supplementary figs. S2 and S4, Supplementary Material 
online) (Zemach et al. 2010; Bonasio et al. 2012). 

The fire ant displayed substantially lower levels of intragenic 
DNA methylation than the honey bee (supplementary fig. SI, 
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Table 1 

DNA Methylation Largely Targets the Same Genes in the Fire Ant and 
the Honey Bee 



Species 


Methylated 


Unmethylated 


Uniquely 




Genes 


Genes 


Methylated in This 








Species (%) 


Fire ant 


2,581 


2,622 


174 (6.7) 


Honey bee 


3,030 


2,173 


623 (20.6) 



Supplementary Material online). There were also minor, but 
statistically significant, differences in coding sequence DNA 
methylation levels among morphs within both species (all 
comparisons between morphs within the fire ant and within 
the honeybee differed significantly, P<0.01, Wilcoxon 
signed-rank tests, n = 5,203 orthologous genes) (Lyko et al. 
2010; Bonasio et al. 2012). 

Comparative analyses using empirical DNA methylation 
data revealed that the localization of DNA methylation was 
highly similar in the fire ant and the honey bee (Spearman's 
rho between species-level DNA methylation levels of ortholo- 
gous coding sequences = 0.71 7, P< 2.2 x 10"^^ n= 5,203; 
supplementary table S4 and figs. S2-S4, Supplementary 
Material online). DNA methylation was targeted to fewer 
genes in fire ants than in honey bees (table 1). However, 
93% of orthologous genes methylated in the fire ant were 
also targeted by DNA methylation in the honey bee (table 1). 
Thus, the genes targeted by DNA methylation in the fire ant 
were largely a subset of those targeted in the honey bee 
(table 1). Furthermore, as in the honey bee and other inverte- 
brates (Hunt et al. 201 0; Sarda et al. 201 2), genes targeted by 
DNA methylation in the fire ant exhibited enrichment of func- 
tional annotations linked to cellular "housekeeping," includ- 
ing an 8-fold enrichment of "translation" annotations 
(supplementary table S5, Supplementary Material online). 

DNA Methylation Is Targeted to Genes Marked by Active 
Histone Modifications in Insects 

The genomic localizations of many histone modifications are 
highly conserved among taxa (Bernstein et al. 2005; Feng and 
Jacobsen 2011; Woo and Li 2012), suggesting that insight 
into histone modification targeting can be gleaned from dis- 
tantly related organisms (Nanty et al. 2011). Moreover, the 
presence of many histone modifications, previously identified 
in diverse eukaryotes, was recently confirmed in the honey 
bee (Dickman et al. 2013). Accordingly, we compared pat- 
terns of DNA methylation in the fire ant and honey bee to 
histone modifications in D. melanogaster (Celniker et al. 
2009). The goal of this analysis was to determine whether 
intragenic DNA methylation was associated with particular 
histone modifications in insect genomes (Nanty et al. 201 1). 

We observed striking positive correlations between DNA 
methylation levels, in the fire ant and honey bee, and several 
histone modifications associated with active transcriptional 



states in D. melanogaster (fig. 1 and table 2) (Filion et al. 
2010; Kharchenko et al. 2011; Nanty et al. 2011). The 
correlation coefficients describing these associations varied 
by exon position according to the spatial localization of each 
histone modification (fig. 1 and supplementary fig. S5, 
Supplementary Material online). In particular, active histone 
modifications that were localized to 5^-gene regions in 
D. melanogaster {e.g., histone H3 lysine 4 di- and trimethyla- 
tion [H3K4me2/me3]) (Kharchenko et al. 2011) were most 
strongly correlated with DNA methylation in 5^-gene regions 
following the translation start site in the fire ant (fig. 1). By 
comparison, active histone modifications associated with tran- 
scriptional elongation (e.g., H3K79me2 and H3K36me3) 
(Kharchenko et al. 2011) exhibited greater homogeneity in 
correlations with DNA methylation across exon positions 
(fig. 1). In addition, the modifications H3K27me2/me3, 
which are associated with domains of polycomb-mediated 
repression and the regulation of tissue-specific gene expres- 
sion (Bell et al. 2010; Kharchenko et al. 2011; Lafos et al. 
2011), were negatively correlated with DNA methylation in 
the fire ant and honey bee (fig. 1). Overall, these correlations, 
which were largely unchanged when controlling for sequence 
length (supplementary table S6, Supplementary Material 
online), suggest remarkable concordance between the pat- 
terning of DNA methylation, when present, and histone mod- 
ifications across evolutionary time (Nanty et al. 201 1). 

DNA Methylation and Active Histone Modifications Are 
Associated with Constitutive Gene Expression 

We investigated the relationship between intragenic DNA 
methylation and gene expression in the fire ant (Ometto 
et al. 2011). We observed a generally positive but bell- 
shaped relationship between levels of DNA methylation and 
overall gene expression levels (fig. 2A), as observed in previous 
studies in diverse eukaryotes (Zemach et al. 2010; Bonasio 
et al. 2012; Jjingo et al. 2012). We further observed a striking 
negative correlation between variability in gene expression, as 
measured by the CV among individuals, and DNA methylation 
(fig. 2B). 

We next investigated the transcriptional correlates of active 
histone modifications in D. melanogaster. Specifically, we ex- 
amined whether particular D. melanogaster histone modifica- 
tions were associated with D. melanogaster gene expression 
measures (Chintapalli et al. 2007) including expression level, 
tissue specificity, and variability of gene expression among 
individuals as measured by CV. We found that active histone 
modifications were positively correlated with overall gene ex- 
pression level and negative correlated with tissue specificity of 
expression in D. melanogaster {f\g. 20, as observed previously 
(Filion etal. 2010; Kharchenko etal. 2011). We also observed 
negative correlations between active histone modifications 
and variability of gene expression among D. melanogaster 
individuals. 
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Fig. 1 . — DNA methylation is targeted to genes marked by active histone modifications. Bars indicate mean DNA methylation level (mCG/CGaii) in the fire 
ant with 95% confidence intervals. The heatmap conveys the strength and direction of Spearman's rank correlation coefficients (rho) between DNA 
methylation levels in the fire ant Solenopsis invicta and enrichment of histone modifications targeted to orthologs in the fly Drosophila melanogaster. Histone 
modifications in green are associated with active transcription and those in purple are associated with repression of transcription (Kharchenko et al. 201 1). 
Exon position is indicated from 5' to 3' relative to the translation start site in each species and "gene" indicates the gene body (introns included but excluding 
untranslated regions). 



Table 2 

Spearman's Rank Correlations (p) between Coding Sequence DNA Methylation Level in the Fire Ant or Honey Bee and Histone Modifications 



Associated with Orthologous DNA Sequences in Drosophila melanogaster 


Drosophila 


Active/Repressive^ 


Genomic Association^ 


Fire Ant 


Honey Bee 


melanogaster 






Methylation (p X)^ 


Methylation (p Xf" 


Histone 










Modification (X) 










H3K4me3 


Active 


Transcription start site-proximal regions 


0.618**** 


0.529**** 


H3K4me2 


Active 


Transcription start site-proximal regions 


0 572**** 


0 505**** 


H3K79me2 


Active 


Transcription start site-proximal regions, transcriptional elongation 
of exonic regions, and intronic regions 


0.527**** 


0.484**** 


H3K9ac 


Active 


Transcription start site-proximal regions 


0.436**** 


0.354**** 


H3K36me3 


Active 


Transcriptional elongation of exonic regions 


0.397**** 


0.392**** 


H4K16ac 


Active 


Chromosome X genes and male dosage compensation 


0.362**** 


0.333**** 


H3K79me1 


Active 


Transcriptional elongation of exonic regions and intronic regions 


0 301 **** 


0 357**** 


H3K4me1 


Active 


Intronic regions 


0.159**** 


0.198**** 


H3K18ac 


Active 


Intronic regions 


0.007^^ 


0.014"^^ 


H3K9me2 


Repressive 


Pericentromeric heterochromatin and heterochromatin-like 
domains 


-0.032* 


-0.01 0''' 


H3K9me3 


Repressive 


Pericentromeric heterochromatin and heterochromatin-like 
domains 


-0.068**** 


-0.025'^' 


H3K9me1 






-0.151**** 


-0.120**** 


H3K36me1 


Active 


Intronic regions 


-0.157**** 


—0 139**** 


H3K27me3 


Repressive 


Polycomb-mediated repression 


-0.316**** 


—0 297**** 


H3K27me2 






-0.432**** 


-0.398**** 



^Kharchenko et al. (2011). 

"^n^ 5,619 orthologs for the honeybee and 4,674 orthologs for the fire ant. 

****P< 0.0001. 

*P<0.05. 

'^^P>0.05. 
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Fig. 2. — Both DNA methylation and active histone modifications are similarly associated with transcriptional regulation. Mean fire ant coding sequence 
DNA methylation level (mCG/CGaii) exhibits (A) a generally positive but bell-shaped relationship with gene expression level, and (B) a negative relationship 
with the CV in gene expression (CV expression). Error bars indicate 95% confidence intervals (n = 3,201 for panels /\ and B). (0 Gene expression measures in 
Drosophila melanogaster are strongly negatively (CV expression and tissue specificity) or positively (expression level) correlated with enrichment of active 
histone modifications in D. melanogaster and DNA methylation in Solenopsis invicta and Apis mellifera (n = 3,892 orthologous genes). All correlations are 
highly significant (P<2.2x 10"^^). 



We next investigated whether DNA methylation levels in 
fire ant and honey bee genes were correlated to gene expres- 
sion nneasuresfor D. me/anogasfer orthologs. Rennarkably, we 
found that nneasures of DNA methylation in both the fire ant 
and honey bee were highly associated with D. melanogaster 
gene expression metrics. Specifically, DNA methylation in the 
fire ant and honeybee was negatively associated with CV and 
tissue specificity but positively associated with expression level, 
in D. melanogaster 2Q. These results support the concor- 
dance between DNA methylation and histone modifications, 
not only spatially (fig. 1) but also functionally (fig. 2Q. 

Promoter Architecture and Intragenic Epigenetic 
Modifications 

Recent studies in mammals and flies revealed that core pro- 
moters of genes fall into distinct classes defined by their RNA 
polymerase II (Pol II) initiation profiles (Lenhard et al. 2012). 
"Broad" core promoters exhibit multiple RNA Pol II initiation 
sites within the promoter and are associated with genes 
exhibiting cellular housekeeping functions and constitutive ex- 
pression (Engstrom et al. 2007; Lenhard et al. 2012). In con- 
trast, "sharp" core promoters exhibit a more focused RNA Pol 



II initiation site and are associated with tissue-specific expres- 
sion (Engstrom et al. 2007; Lenhard et al. 2012). 

To explore potential associations between core promoter 
type and intragenic DNA methylation, we compared coding 
sequence DNA methylation levels in the fire ant and honey 
bee to promoter classifications of D. melanogaster orthologs 
(Hoskins et al. 2011). We found that fire ant and honey bee 
orthologs of D. melanogaster genes with broad promoters 
showed considerably more intragenic methylation than ortho- 
logs classified as having sharp promoters (fig. 3 and supple- 
mentary fig. S6, Supplementary Material online). 

Discussion 

Intragenic DNA Methylation and Active Histone 
Modifications 

The targeting of DNA methylation appears to be highly con- 
served across invertebrate taxa that possess functional DNA 
methylation systems (table 1) (Sarda et al. 2012). This target- 
ing also appears to be strongly associated with the presence of 
several active histone modifications in insects (fig. 1 and 
table 2). Moreover, DNA methylation is preferentially targeted 
to 5^-regions within insect genes (fig. 1) (Zemach et al. 2010; 
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Fig. 3. — DNA methylation patterns differ based on promoter archi- 
tectural classification. (A) Mean DNA methylation level in fire ant coding 
sequences differ according to sharp (n = 530) and broad {n = 2,91 9) core 
promoter architecture classifications of Drosophila mela nogaster orthologs 
(Mann-Whitney Latest, P< 2.2 x 10"^^). (B) Density graph comparison of 
fire ant DNA methylation distributions for each promoter class-based gene 
category, illustrating high relative representation of methylated genes as- 
sociated with broad core promoters. 



Bonaslo et al. 201 2). H3K4me2/me3, which is associated with 
5^-gene regions in D. me/anogasfer (Kharchenko et al. 201 1), 
exhibits the strongest correlations with DNA nnethylation in 
our analyses (fig. 1). 

The associations between DNA methylation and H3K4me2/ 
me3 are particularly rennarkable because H3K4 methylation is 
negatively correlated with intragenic DNA methylation in 
mammals (supplementary fig. S7, Supplementary Material 
online) (Ooi et al. 2007; Meissner et al. 2008; Maunakea 
et al. 2010). One possible explanation for this difference is 
functional divergence of DNA methyltransferase enzymes 
between insects and mammals (Ooi et al. 2007), facilitating 
differential preference to methylated versus unmethylated 
H3K4 in the two taxa. Regardless, in contrast to mammals, 
intragenic DNA methylation appears to be consistently posi- 
tively associated with a wide array of active histone modifica- 
tions in insect genomes (fig. 1, supplementary fig. S7, 
Supplementary Material online; see supplementary discussion. 
Supplementary Material online). 

Intragenic DNA Methylation and Gene Regulation 

Intragenic DNA methylation is targeted to genes with low 
variability of expression among fire ant individuals (fig. 2B). 
DNA methylation is also associated with genes displaying 
ubiquitous expression among tissues in the honey bee (Foret 
et al. 2009) and morphs in ants (Bonasio et al. 201 2), as well as 
with orthologs that exhibit ubiquitous expression among 
D. me/anogasfer tissues (fig. 2Q. These results, coupled with 
the broad core promoter architecture of orthologs of methyl- 
ated genes (fig. 3), provide strong evidence that intragenic 
DNA methylation is preferentially targeted to constitutively 
expressed genes in insects. 



DNA methylation is known to influence transcriptional reg- 
ulation in a number of ways. For example, DNA methylation 
may directly alter transcription factor binding (Bird 2002; 
Shukia et al. 201 1 ; Wang et al. 201 2). In addition, DNA meth- 
ylation may increase the compaction and rigidity of nucleo- 
some-DNA complexes (Choy et al. 2010; Lee JY and Lee T-H 

2011) , thereby indirectly altering DNA accessibility to tran- 
scription factors (Li and Widom 2004; Bintu et al. 2012). In 
this manner, intragenic DNA methylation may shield cryptic 
binding sites within gene bodies, thus preventing spurious 
transcription initiation in constitutively expressed genes 
(Zilberman et al. 2007; Maunakea et al. 2010). Intragenic 
DNA methylation is also preferentially targeted to exons (sup- 
plementary figs. S2 and S3, Supplementary Material online), 
where RNA Pol II is most prevalent (Yin et al. 201 1) and nu- 
cleosomes are preferentially positioned in diverse eukaryotes 
(Schwartz et al. 2009). This underscores the possibility that 
intragenic DNA methylation acts to enhance or modify tran- 
script integrity and nucleosome positioning in insects 
(Chodavarapu etal. 2010). 

The relationships between DNA methylation, histone mod- 
ifications, nucleosome positioning, and histone variants 
(Coleman-Derr and Zilberman 201 2) suggest that DNA meth- 
ylation is complementary, and potentially redundant, to other 
forms of epigenetic information in eukaryotes. We speculate 
that the multilayered nature of the eukaryotic epigenome may 
help to explain how the loss of DNA methylation has been 
tolerated in some insect lineages (Zemach et al. 201 0; Glastad 
etal. 2011). 

Histone Modifications, DNA Methylation Patterning, and 
Gene Expression Profiles Define Constitutive and Variable 
Transcription States in Insect Genomes 

The association between the patterning of DNA methylation 
and active histone modifications (fig. 1), along with their sim- 
ilar regulatory relationships (figs. 2 and 3), suggests that genes 
can generally be partitioned into two highly conserved epige- 
nomic states (fig. 4). These states can be described as "con- 
stitutive" and "variable" (fig. 4). Constitutive and variable 
states differ in transcriptional activity (Filion et al. 2010; 
Kharchenko et al. 2011; Hunt, Ometto, et al. 2013), pro- 
moter architecture (Hoskins et al. 2011; Lenhard et al. 

2012) , epigenetic modifications (Foret et al. 2009; Filion 
et al. 2010; Kharchenko et al. 2011), physical domains of 
the genome (Sexton et al. 2012), gene function (Filion et al. 
2010; Sarda et al. 2012), and selective constraint (Sarda et al. 
201 2; Hunt, Ometto, et al. 201 3). The fact that we were able 
to infer constitutive and variable states from analyses of DNA 
methylation data from Hymenoptera (5. invicta and A. melli- 
fera) and diverse data sets from Diptera (D. mela nogaster), 
two insect orders that diverged approximately 350 Ma 
(Wiegmann et al. 2009), suggests that such regulatory 
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Fig. 4. — Histone modifications, DNA methylation patteming, and gene expression profiles define constitutively and variably expressed genes in insects. 
A heatmap of coding sequence DNA methylation levels, histone modification enrichment, and gene expression measures produced by the Ward hierarchical 
clustering method with dendrograms proportional to distance. Two clusters emerge (n = 3,892 orthologous genes). 



properties nnay be highly conserved over evolutionary tinne 
(Engstrom et al. 2007; Hunt, Glastad, et al. forthcoming). 



Conclusions 

intragenic DNA nnethylation exhibits limited variation among 
insect taxa with functional DNA methylation systems and is 
strongly associated with multiple histone modifications tar- 
geted to orthologous D. melanogaster loci. Many regulatory 
properties of DNA methylation are potentially driven, at least 
in part, by the association of DNA methylation with conserved 
epigenetic and regulatory genomic domains, illustrating that 
DNA methylation is a single component of a conserved, 
complex, multilayered epigenome. Overall, our results support 
the view that the regulatory roles of intragenic DNA methyl- 
ation cannot be fully understood without considering the tight 
integration of DNA methylation with histone modifications, 
nucleosome positioning, and RNA Pol II kinetics. 
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Supplementary materials and methods, discussion, figures 
S1-S7, and tables S1-S6 are available at Genome Biology 
and Evolution online (http://www.gbe.oxfordjournals.org). 
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